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ABSTRACT 


“21, INTRODUCTION 


The cardiovascular system, which refers to the heart, blood vessels, and the blood is the basis of the vivacity of 
virtually all animated life. The heart (cardio) ensures the pumping of blood and the blood vessels (vasculature) are 
the appurtenances that passively maintain the integrity of flow. The flow of blood in the vessel is driven by the 
pressure difference between any two arbitrary points along the length of the vessel. While the latter is true, the 
pressure difference across an organ is described by the variation between arterial pressure (P,4) and venous pressure 
(Py). The circulatory system is, by default, marked by laminar blood flow in which concentric layers of blood flow 
parallel down the length of a blood vessel. However, in more or less pathological conditions, the laminar flow 
regime may be upset by turbulence in which the critical Reynolds number (Re) is exceeded. Turbulence is, to all 
intent and purposes, chaos. Both the vasculature and the cardiac chambers may experience distressed regions being 
implicated or exacerbated by turbulent flow. In all, pressure fluctuations and their incidences are produced by an 


ensemble of fluid-structure interaction-based factors. 


In essence, the pulsating cardio-vascular structures, distressed areas of turbulent flow, mixing of flows of varying 
temperatures, and their semblances produce pressure fluctuations that are transmitted through the fluid as sound. In 
both the vasculature and the cardiac chambers, this sound due to turbulent flow is the progression of compressions 
and rarefactions propagated away from the source environs. In the former, such sounds are called bruits such as can 
be auscultated in the carotid, renal, and visceral arteries among others. Stenosis is a culprit in diverse abnormal flow 
conditions in the vessel. Seo and Mittal [1] suggested that the chief source of arterial bruits is perturbation initiated 
by vortex flow in the near post-stenotic region. Lees and Dewey [2] proposed phono-angiography as a non-invasive 
diagnostic technique. In Fredberg [3] the fluctuating mechanical stress on the arterial-wall downstream of stenosis 
by turbulent blood flow was investigated. The work found that the distance downstream of the stenosis at which the 


mean-square fluctuating pressure attains its maximum intensity was dependent upon jet velocity and fluid viscosity 
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only. In an elegant study Sack et al. [4] advanced a “sound boundary” curve which creates the combinations of 
minimum flow Reynolds numbers and per cent stenosis requisite for the onset of vascular murmurs. As reported in 
[4], Bruns [5] argued that periodic vortex shedding created at the site of vascular obstruction in the form of vortex 
rings is the source of vascular murmurs- a claim corroborated by Fruehan [6]. However, in the case of an orifice, 
vortex shedding speedily acquiesces to the instabilities of the jet and eventual turbulence at very low Reynolds 


numbers [7]. 


Heart murmurs are defined concerning the location in the thoracic cage, pitch, volume, and phase of the cardiac 
cycle [8]. They are often explained in association with the heart valves. Each valve opens or closes in accord with 
differential blood pressure on each side [9, 10] and allows flow in one direction through the chambers of the heart. 
Vascular bruits or heart murmurs, whichever sound may be blamed on the stenosis or partial occlusion of the blood 
vessels or the valve. It should be noted that functional murmurs (also called physiologic murmurs (?) (see 
Napodano [11]), which emanate from sources away from the heart can ensue in the absence of valvular 
pathology. Sabbah and Stein [12] detailed studies indicating that turbulent flow takes place above the aortic valve 
under some circumstances in normal humans much as it typifies the incidence of aortic valvular disease. Becsek et 
al. [13] depicted free-stream instabilities of the core jet in addition to instabilities related to fluid-structure 
interaction between blood flow, as some hydrodynamic instability, leading to turbulent blood flow in the wake of 
the aortic valve. In the main, the critical issue regarding bruits and murmurs is the sound source. Acoustic 
auscultation for cardiac sounds is largely done in the regions where the cardiac sounds are best heard- aortic, 
pulmonic, tricuspid, and mitral [14]. However, these regions do not relate to the location of heart valves. 
Since Doppler techniques reveal insignificant regurgitation across atrioventricular valves that do not produce 
audible murmurs [12], Tavel [16] proposed a well-informed identification and classification of audible systolic 


murmurs, which would facilitate the detection of their possible source and mechanism. 


“<2, FLOW-INDUCED VASCULAR SOUND 


Sound arises from the transmission from pressure fluctuations induced by the pulsating cardio-vascular structures, 
distressed areas of turbulent flow, mixing of flows of varying temperatures, and their semblances. In both the 
vascular- ture and the cardiac chambers, this sound is the progression of compressions and rarefactions propagated 
away from the source environs. Vascular sounds may be normal, innocent findings (e.g. a venous hum in a child) or 


may point to underlying pathology (e.g. a carotid artery bruit caused by atherosclerotic stenosis in an adult). 


2.1. Stenosis in the vascular compartment 


Figure 1. Flow through a intact and stenosed vascular compartment. D is the diameter of the intact vessel, while d 
is the diameter of the stenosed vessel 
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In the course of systole, blood proximate to a stenosis in a laminar flow experiences rapid convective acceleration 
as it passes from the intact compartment of the artery through the converging section of the stenosis. Owing to the 
inability to overcome the adverse pressure gradient, the fluid passing through the diverging section of the stenosis 
separates from the walls. This occurs at a sufficiently high Reynolds number. A shear layer is produced at the 
boundary between the separated jet, which is at a high velocity, and the relatively slower-moving fluid in the 
recirculating separation zone, which is amenable to fluid-dynamical instabilities. The shear layer created affords a 
source from which such instabilities educe energy from the mean flow [17, 18]. The energy extraction process 
progresses at a considerably rapid rate that the instabilities culminate in a fully turbulent motion before the end of 
systole if the jet Reynolds number is sufficiently high. Once the jet fills the artery, the decay of turbulence results as 
it (turbulence) will no longer sustain itself by extraction of energy from the mean flow because the intact arterial 
Reynolds number is characteristically below the critical Reynolds (Reégit ~ 2000) number required to attain 


sustained turbulent flow in a straight arterial configuration. 


In describing the vascular flow we assume that the blood is a Newtonian (see Pedley [19]), there is no vessel wall 
deformation in the flow field, the acoustic waves dominate the shear waves generated in the tissue, and the viscous 
dissipation of the acoustic wave is negligible. In the event of stenosis the equation of constriction in a stenosed 


arterial compartment is of the form (see [1, 20]) 


() 


Y=Voue e 1+ o0s{ 20 == » -H<(x-x,)<H 
2 H 


where / is the size of constriction, xp is the centre of the stenosis, and His the height of the vascular tissue (assumed 
homogeneous). However, the axisymmetric mathematical problem may be given by the two-dimensional radial, r, 


and axial, z, co-ordinates, as detailed in Borisyuk [21] such that 


r=r_ (for the intact vessel) (2) 


0 


r=Ar, (for the narrowed vessel) 


where 7” = A/Ag encodes the ratio of the minimum cross-sectional area, A, to the open (non-stenosed) lumen area, 
Ao, of the vessel (0 <A’ < 1). The equivalent variations in the mean flow characteristics are considered from the 


situation of the same flow in the narrowed and intact arteries. Thus, 

uA = UAg, (3) 
with the mean flow velocity, u, in the stenosis given by 

u = U (Ap/A) = U/ 2. (4) 
Arterial turbulence is assumed in a flow with Reynolds number Re, 9 = U (2r9)/v in which condition Re,» > Re¢rit 
holds [22, 23]. For flow in a pipe of diameter D, experimentation indicates that for “fully developed” flow, the 


critical Reynolds number is about Re .,i, = 2300 [24]. However, the development of turbulence in the (large) blood 


vessel is taken as 2000 [25, 26, 17]. 


2.2. Acoustic analogy and theory 


It is customary that fluid dynamical phenomena obey the laws of conservation of mass and momentum. The generic 


equation relating the former reads: 
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where p and u encode the fluid density and velocity, respectively. The equation of conservation of momentum 


reads: 


20 ¢u.vyn = (Pe 4 Py Fe” - 
Ot on a) x Ox 


i i 


p 


where P;;is the stress tensor. It is of note, from the above two equations, that the use of fluid density is emphasized. 
Density fluctuation is the phenomenon of essence in sounds produced by unsteady liquid flow, whilst pressure 
fluctuation applies more to unsteady combustion (see Stanko [27]). The celebrated “Lighthill analogy” which is an 
inhomogeneous wave equation is of the form 


op' r op! ze a (pu, -9,) 2 (p'-cp') ‘ (7) 
or * Ox? ax,Ox, Ox? 


In the above: c, is the sonic speed in the reference quiescent state (9, po) of the fluid adjoining a listener, p and o 
encode the pressure and viscous parts of the stress tensor, the primed quantities, p’ = p—po and p’ = p—po are 
deviations from this reference states of density and pressure respectively. The generic form of Lighthill's equation 


relating density to turbulence stress tensor, T;;, for the acoustic field is 


2 2 8 
a ~ a = T,, 
ot oxox," 
where 
T, = pu, —0, +(p—c,p)6, ° (9) 


and 6 is the Kronecker delta. Lighthill’s derivation is characterized by an unsteady flow in an unbounded domain, 
ij 


but it has the beauty of a closed-form solution. This solution is [28] 


(10) 
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where V encodes the far-field source region, where |x — y| is the distance between the acoustic observation point, 
x, and the source field/point, y, n is the outward normal from the fluid. This quadrupole field may be written in a 


more simplified form as 


= es 2 : : i 
1, @wW-9)1 SF fy, PW ay ae 
Ix-y| c, Ot” c 


0 


As glorious as Lighthill’s derivation (analogy) is, its shortcoming is the absence of flow boundaries. Each 
component of the vasculature (arteries, arterioles, veins, venules) consists of flow boundaries, all with 
characteristics which influence flow patterns. Fortunately, Curle [29] bridged this gap by showing that boundary 
terms could largely bolster the mass and momentum of the flow. It is pertinent to appreciate the creation of better 
monopole and dipole sources of profoundly greater acoustic strength in Curle’s creation; Lighthill’s was much of 


volume quadruples. In the presence of rigid and stationary surfaces in the flow, the solution to (8) reads [30]: 
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where G = G(x, t;y,T) is Green’s function, and n; are components of the outward unit normal of the surface S, 
pointing into the fluid. Using free-space Green’s function, Curle’s solution to the Lighthill equation in the acoustic 
far field takes the form: 


1 rr, r 1 ro r (13) 
x,t) © ‘i _T') y ¢ d’y+ stn t—— |d’ 
pt) 4nc: I Lp £) = Arc’ Jos Ot snl £) 3 
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where r = |x — y|,7; = xX; — yj. Curle’s surface integral predicts a compact dipole when the solid body is 


negligibly small. 
2.3. Vascular sound field from a source point 


In the course of reading this paper, it is instructive to emphasise that whereas bruits (murmurs) associated with 
turbulent flows are largely pathologic, a good many cases are common in healthy individuals. The epigastric 
murmur-a systolic medium-to-low-pitched murmur over the epigastrium [31], and cervical bruits and hums which 
may arise from neck arteries or veins [32] are a few cases of murmurs that may present innocuous findings or 
indicate underlying pathology. From the foregoing, a non-stenosed atrial compartment subject to turbulence is, in 
all possibility, amenable to bruit. Here we consider a single point source in the arterial compartment as a source 
whose dimensions are much less than the wavelength of sound generated. Details regarding turbulent fields in short 
vessels are sparse, and therefore an approximation of a random point source may suffice. Suppose x9 € (d), /2) is a 


source point in the arterial segment. The defining model of the source point is of the form 
Pru (a, X,) = Pu (o)1 O(x = x,) (14) 


where Pj, is the turbulent pressure, w the angle frequency, ] the effective length of the source, and x, is the 
location of an equivalent source. Here the 6 function has units of L’' (see Succi [33]). Applying the spatial Fourier 


expansion of the left-hand side of (14) gives 


1 + : 
Pan lOrK)= Fe | APs (Oa Jem (15) 
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where «K encodes the wave number. (We note that in the arterial segment considered, the lower limit of integration, 
l;, may be considered as the pipe origin at the inception of turbulence; in that case, it may be numerically zero). It 
will be reasonable to assume that the turbulent sources are contained entirely inside the pipe and are sufficiently far 
from the ends. Therefore, p,,,, and the convective derivative dpj,p/dxo could be set to zero at the inlet and exit. As 


well, assume that the observation point is situated outside the source region. The acoustic field of the point source is 


. 16 
panty =Paa@rbdiad | 1 peng Lp ginnrngrnn | (16) 
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where M encodes the Mach number, R;, pressure reflection coefficient at pipe exit, /,. The radiated field is 


(16a) 


< G,(x| X,)s 


0 


P(@,X) = iN Pi, (Q, “fr —iM 
where G is Green’s function which gives the linear acoustic response of a pipe. The pressure field induced by a 


harmonic source of unit strength located at x = xo is Green's function G,(x|xo ). By getting the magnitude | p(a,l, yf 


in (16a) the averaging over all source positions removes all dependence on x,. Thus, 


(17) 


2) Le > |p ,(@)| wre? 1 1 2 
1) )=—f' ds 1), = . 7 =|R, 
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In a similar mode, the average of the acoustic field of a distributed source over N measurements of incident pressure 


is shown [33] to be 


2 12 2 
(Pot) ) =X iP@.1, (18) 
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£3, VASCULAR ACOUSTIC POWER 


The hydrodynamic near-field is the region contiguous to the vibrating surface of the source. Vascular bruits 
produce vibration. Interestingly, in the near-field region, the fluid motion is not directly related to sound 
propagation and measurements of the acoustic pressure amplitude do not indicate the sound power radiated to the 
far field by the source [34]. Radiated sound power is predictable from an appropriate number of pressure 
measurements made in the geometric (Fresnel) near field. The far field, also known as the radiation-zone field, 
carries a somewhat uniform wave pattern. The radiation zone is essential since far fields largely fall off in 
amplitude by 1/r. Thus, the total energy per unit area at a distance r is proportional to 1/r’. As the area of the sphere 
is proportional to r’, the total energy passing through the sphere is constant. The far-field energy therefore radiates. 
Far-field strength decreases with increasing distance from the source. This results from the so-called inverse-square 
law for the radiated power intensity of radiation. A source (point) is considered to be in the far field if the distance 
between the source and a measuring device is beyond a few times the maximum acoustic wavelength and the 
device’s dimensions are small compared with the minimum wavelength, with basically planar inward bound 


wavefront at each device. The stethoscope is the auscultation (measuring) device in this case. 
3.1. Power spectrum 


We now consider the fields in line with vascular sound radiation. The acoustic power spectrum is considered the 
basis for studying noise fields produced by physiologic flows [35, 36], 37]. Fredberg [17] underscored this by 
having investigated the spectral density of pressure as one of the main descriptors of the fluctuating turbulent 
pressure at the vascular wall. The power spectral density (PSD) of the pressure may be gotten from the fluctuating 
Poisson equation (see Lee ef al. [38], Panton and Linebarger [39]). For an incompressible flow, the resulting 


Poisson equation for the fluctuating pressure reads [40] 
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where summation over repeated indices is implied, the lower case quantities with ’ (prime) indicate fluctuations, 
and the upper case quantities with bars indicate mean quantities. Over bars in general indicate time-averaging. The 
first right member of (19) above, a linear source term, encodes the mean-shear-turbulence interaction; the second 


member, a non-linear source term, encodes the turbulence-turbulence interaction. 


In the light of the Lighthill acoustic analogy, consider the space-derivative form 


a (20) 
Pay (Xt) = —|— y be 
4a 7 


For a far-field location of the observer point with |x| >> [y| the integrand in (20) may be approximated by the second 


time derivative taken at the retarded time (see Sarkar and Hussaini [40]). Thus, 


1: xe postal) y (21) 
4nc, x r ; 


Pay (% t) = 


where 7, (y,t—r/c,)indicates 9°(pyy)/ ar evaluated at retarded time ¢ - r/cp. In using a two-point correlation 
T (0) T ( r) from the (Direct Numerical Simulation) DNS acoustic power may be obtained. 


In line with the above, suppose (x, 0, 0) is an observation point on an arterial segment. The acoustic power P,(f) per 


unit mass at the observer point given by 


P(t) =——-||T_(&,,t,)T_(x,,t, dx dx (22) 
( ) AnVc ff me A a) mi B jd A B 


where V encodes the source volume, and X, and f, encode the position of source point A and time of sound radiation 


at point A to reach the observer at time ¢; similarly, X, and t, apply to source point B. If, in line with Proudman [41], 


T (x, tT (x,t) is assumed constant over the maximum time difference fy - tg thereby disregarding the retarded 


time effects when integrating the volume with respect to source point B one finds equation (22) in the form 


1 - = 
P.()=— JT. tT Apt ax, ° oy 
Arc; , 
Proudman [41] further considered the case where t, is constant in the course of the integration over source point A. 
Thus, 


1 
Atc 


P(t)=——JT (x, )T (x, )dx,” — 


Note that the two-point spatial autocorrelation in (24) is a function of r = |x4-x,| only, owing to the turbulence 


isotropy. From da Cunha Lima ef al.[43], the power spectral density may be employed in examining possible 
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alterations in a frequency distribution as the observation point moves downstream from obstruction in a vessel. 


Thus, the power spectral density of the turbulence at a given location r may be given by 
P(r, @) =v(r,@) v(r, @), (25) 
in which 
v(r,@) = [0 v,ar,pe" dt. (26) 
where v(r, @ = 2zf) is the Fourier transform of the x component of the velocity fluctuation at location r, 


corresponding to w = 2af, f encodes the frequency. The key interest is in the power spectra corresponding to the 


x-direction contribu- tion to the energy. 


With the velocities specified at N discrete values of time, at time interval At, the discrete Fourier transform is 


calculated at frequencies 
Vv, = (r, f, = n/ NA) => Sv(r, t, er . 5 
k=0 


At a given position r inside a frequency interval J the power of the fluctuation is the integral 


; 28 

AP(r) = J, P(r, f df ~A*Y'\v, 0) oo 
where M\v,(r)/° = P(r, @ = 2zf) is the discrete Fourier transform of the power spectral density (see [43] ). Note that 
the quantity inside the summation in (28) has the dimension of energy. The acoustic power that issues from a 


volume V of fluid reads [40] 


2 2 
p(t) =P Oe (29) 


PK 
and the acoustic power per unit mass of turbulent fluid is 


p, (t) (30) 
pV 


0 


PO= 


In (29) Dims encodes the root-mean-square pressure. We have 


a ec = 31 
Pirie Pa) = (> if pata oe 
where T encodes a suitable averaging time. The spectral density pressure is in the form 
Prortg, = Jo ECS Ff - (32) 


For the radiated sound in the far field the intensity /,(x, t) is proportional to the square of the fluctuating pressure 


Purv. This intensity is defined by 


| 33 
1,(x,t) = (Bue P) ° 
PG 
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where pp and p are the fluid density and pressure, co is the ambient sound speed. The intensity of radiated sound at 


any point x, and at time ¢ can be determined by [44, 45] 


1 (34) 


fou 
1.(x,t) = ——— [|| dy [Jar -—P, (y,r,7). 
=e) oy lllavillor P.o.80 


In the above y is a position, r is the distance vector between two internal points y, and yz (i.e. r = |y4-ya|) inside the 
integration volume and T is the retarded time, x =(x|. The space-time correlation function P,,., can be written in the 


form 


P.=AUw'(r/ DOT) (35) 


where u, / and [ are, respectively, reference values of the turbulent velocity, the turbulent length scale, and 
turbulent frequency, (1/1) is the non-dimensional space correlation function, Ap = constant is a scale coefficient, 


and @([‘r) is the non-dimensional retarded time correlation function. The modes for the space and temporal 


correlation functions adopted in Lilley [46] enhanced Proudman’s estimate of the acoustic power per unit of 


volume in the time domain in the form 
P, =A,peM; (36) 


where M, encodes the turbulent Mach number. The analysis based on two primary descriptors of the fluctuating 
turbulent pressure at the vascular wall: the root-mean-square pressure and the spectral density of pressure, which 
were investigated [17] indicates that the pitch of vascular murmurs is proportional to the mean systolic velocity in 


the intact portion of the artery and inversely proportional to the stenosis diameter. 
"4, DISCUSSION 


This work looked at vascular murmurs from a generalized theoretical standpoint, noting that the vasculature 
consists of several vascular tissues. It is to be noted that physiological murmurs due to blood flow across normal 


valves exist and may be amplified by both isotonic and submaximal isometric (handgrip) exercise (see [48]). 


Figure 1 shows the flow through a stenosed vascular compartment. The cylinder with open ends A (inlet) and B 
(exit) is the non-occluded arterial lumen in the stenosed region, hypothesised as a mini-arterial compartment. The 
flow upstream (towards the stenotic region) is laminar. The inception of stenosis, is marked by a high fluid velocity 
through the inlet, A. As the velocity increases, shear flows go through a sudden transition from laminar to turbulent 
motion and the inception of turbulence radically changes transport efficacy and mixing properties. A jet with 
Reynold's number, Re;= u;d/v (where v is the kinematic viscosity of the blood) ensues at the orifice, B. Turbulent 
pressure field at the arterial wall transmit oscillatory stress to the surface of the skin, resulting in substantial 


alterations of the spectral density of the fluctuating stress (see Fredberg [47]). 


The energy being extracted by turbulence from the mean flow persists as the jet develops to fill the artery. 
Subsequently, the turbulence no longer abides because the intact arterial Reynolds number (Re = Red/D)) is 
characteristically below the critical Reynolds number required to achieve continued turbulent flow the pipe 


(assumed straight, with Re ~ 2000). 
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From the foregoing discussion: 


(a) One would expect some whooshing sound at the inlet, A (Figure 1) due to the sudden reduction in the arterial 


diameter (D > d), which leads to distortion of the flow. 


(b) The observation location should be predicted within the near field of the sound source. It is pertinent to admit 
that auscultation of bruits sometimes can be felt in a rather far-field, as evident in the aorto-iliac disease that give 


rise to the femoral and mid/lower abdominal bruits. 


(c) A retrograde wave would be expected with the inlet as the source point due to the very high fluid velocity there; 


however, this wave evanesces as it is enveloped by the advancing superior anterograde wave. 


(d) Given (c) above, it is sanguine to posit that incipient vascular murmur pervades from a distance well before the 
source point. The inception of this incipient murmur shall be of clinical importance; it would serve as a guide to 
hazard a guess as to the location of a bruit. But how would this inception be known? The ability to specify the point 


at which the antero-retrograde waves coalesce that leads to evanescence of the retrograde wave furnishes the link. 


(e) The inherent wall shear stress which fluctuates with time and changes course (scouring the arterial wall) has 
spatial gradients which is capable of detaching one endothelial cell from another (local stretching and 


compressing), along with generating structural fatigue within the arterial wall. 


In conclusion, vascular bruits are a source of worries as they are often implicated in cardiovascular events. Postural 
manoeuvres may account for physiological bruits but the notable pathology underlying bruit is vascular stenosis 
which brings to bear on the flow characteristics (the creation of turbulence). Vascular sounds are generated locally 
in peripheral blood vessels and require near-field auscultation. However, retrograde waves that travel at a distance 


away from the main acoustic source create sound on impact with the advancing anterograde wave; this may call for 


a far-field auscultation. 
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